Brown dwarfs forming in discs: where to look for them? 



Dimitris Stamatellos I a and Anthony Whitworth 1 

School of Physics & Astronomy Cardiff University, 5 The Parade, Cardiff, CF24 3AA, UK 



Abstract. A large fraction of the observed brown dwarfs may form by gravitational fragmentation of unstable 
discs. This model reproduces the brown dwarf desert, and provides an explanation the existence of planetary- 
mass objects and for the binary properties of low-mass objects. We have performed an ensemble of radiative 
hydrodynamic simulations and determined the statistical properties of the low-mass objects produced by gravita- 
tional fragmentation of discs. We suggest that there is a population of brown dwarfs loosely bound on wide orbits 
(100 - 5000 AU) around Sun-like stars that surveys of brown dwarf companions should target. Our simulations 
also indicate that planetary-mass companions to Sun-like stars are unlikely to form by disc fragmentation. 



1 Introduction 

The formation of brown dwarfs (BDs) is not well under- 
stood ([35 1). Low-mass objects are difficult to form by grav- 
itational fragmentation of unstable gas, as for masses in the 
BD regime (J 80 Mj, where Mj is the mass of Jupiter) a 
high density is required for the gas to be Jeans unstable. 
Thus, if BDs form according to the standard model of star 
formation, i.e. the collapse of a prestellar core, then the 
pre-BD core must have a density of Z 10~ 16 g crrT 3 . 

E4l and lfT4ll suggest that these high density cores can 
be formed by colliding flows in a turbulent magnetic medium. 
However, this model requires a large amount of turbulence, 
and has difficulty explaining the binary properties of BDs. 
Additionally, the large number of brown-dwarf mass cores 
that the theory predicts have not been observed. 

Another way to reach the high densities required for the 
formation of BDs is in gravitationally unstable discs ([34|, 
1 28 1). These discs form around newly born stars and grow 
quickly in mass by accreting material from the infalling 
envelope. They become unstable if the mass accreted onto 
them cannot efficiently redistribute its angular momentum 
outwards in order to accrete onto the central star (|2|). 
This mechanism reproduces many of the observed statis- 
tical properties of low-mass objects that are not satisfacto- 
rily explained by other theories, in particular the BD desert, 
the statistics of low-mass binary systems, and the forma- 
tion of free-floating planets (ED), ED, El)- 

BDs are also thought to form as ejected embryos from 
star forming regions, i.e. as by-product of the star forma- 
tion process ( I1251D . In this model BDs form the same way 
as low-mass stars, i.e. in collapsing molecular cores, but 
shortly after their formation they are ejected from their 
parental core and stop accreting any further material. Hence, 
they do not realise their potential to become hydrogen- 
burning stars. 



a e-mail: D.Stamatellos@astro.cf.ac.uk 



The above BD formation mechanisms are not exclusive 
of each other and can work in conjunction in star form- 
ing regions. For example, turbulent fragmentation can pro- 
duce cores which collapse to form discs, the discs may then 
fragment to produce BDs, and the BDs may avoid accret- 
ing additional mass by being ejected. 

In this paper we focus on the mechanism of BD forma- 
tion by fragmentation of gravitationally unstable discs. We 
discuss the properties of the BDs produced with this mech- 
anism and the predictions of the model that can be tested 
by observations. 



2 Brown dwarf formation by fragmentation 
of gravitationally unstable discs 

Discs can fragment if (i) they are massive enough that grav- 
ity overcomes the thermal and the local centrifugal support 
([33 1), and (ii) they can cool fast enough that the energy 
provided by the collapse of a proto-fragment is radiated 
away so that the growth of fragment continues ( ifTTIl . |26|). 
Analytical and numerical studies have shown that the cool- 
ing time must be on the order of the dynamical time which 
happens to be similar to the orbital time. 

We have performed an ensemble of radiative hydrody- 
namical simulation of large unstable discs. The evolution 
of the disc is initially followed using Smoothed Particle 
Hydrodynamics (SPH), until ~ 70% of the disc mass has 
been accreted, either onto the stars condensing out of the 
disc, or onto the central primary star; this typically takes 
10 to 20 kyr. Then the residual gas is ignored and the long 
term dynamical evolution of the resulting ensemble of stars 
is followed up to 200 kyr, using an N-body code. 

The energy equation and associated radiative transfer 
are treated with the method of [27 1. This method takes into 
account compressive heating or expansive cooling, viscous 
heating, radiative heating by the background, and radiative 
cooling. The method performs well, in the optically-thin, 
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Fig. 1. SPH radiative hydrodynamic simulation of the evolution of a 0.7 M Q disc around a 0.7 M Q star. The disc is represented using 
2 x 10 6 SPH particles. 3 snapshots are presented. The disc is gravitationally unstable and can cool efficiently, hence it quickly fragments 
to form low-mass H-burning stars, BDs, and planetary-mass objects. These objects form at radii from ~ 100 to ~ 300 AU but due to 
mutual interactions a large fraction of them escape from the disc. 



intermediate and optically-thick regimes, and has been ex- 
tensively tested fll27ll ). The gas in the disc is assumed to be 
a mixture of hydrogen and helium. We use an equation of 
state (| 4], [21 1) that accounts for the rotational and vibra- 
tional degrees of freedom of molecular hydrogen, and for 
the different chemical states of hydrogen and helium. We 
assume that ortho- and para-hydrogen are in equilibrium. 
For the dust and gas opacity we use the parameterization 
by 1 3 1, K(p, T) = Ko p" T b , where kq, a, b are constants that 
depend on the species and the physical processes contribut- 
ing to the opacity at each p and T. The opacity changes 
due to ice mantle evaporation and the sublimation of dust 
are taken into account, along with the opacity contributions 
from molecules and H ions. 

We assume a star-disc system in which the central star 
has initial mass M\ = 0.7 M©. Initially the disc has mass 
M D = 0.7 Mq, inner radius R m = 40 AU, outer radius 
R OVT - 400 AU, surface density E(R) oc /?~ 7//4 , tempera- 
ture T(R) oc R~ l l 2 , and hence approximately uniform ini- 
tial Toomre parameter Q ~ 0.9. Thus, the disc is at the 
outset marginally gravitationally unstable. The radiation of 
the central star is taken into account. 

Three snapshots from a typical simulation are shown in 
Fig. Q] The gravitational instabilities develop quickly and 
the disc fragments within a few thousand years. Typically 
5-10 objects form in each disc, most of them BDs. The 
final stage of the system is determined by the subsequent 
accretion of material from the disc and mutual interactions 
that these objects undergo in the disc. 



3 The properties of objects formed by disc 
fragmentation 
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Fig. 2. The mass spectrum of the objects produced by disc frag- 
mentation. Most of these are BDs (70%); the rest are low-mass 
stars. The vertical dotted lines correspond to the D-burning limit 
(~ 13Mj) and the H-burning limit (~ 80Mj). The red dashed line 
refers to a low-mass IMF with AN /AM oc M~°- 6 (see text). 



3.1 The mass spectrum of the objects produced by 
disc fragmentation 

Most of the objects (~ 70%) formed in the discs are BDs, 
including a few planetary-mass BDs. The rest are low- 
mass hydrogen-burning stars. The typical mass of an object 
produced is ~ 20 - 30 Mj. The shape of the initial mass 
function (Fig. |2]i is similar to what observations suggest, 
i.e. it is consistent with AN/AM oc M~ a , where a * 0.6 
(e.g. in Pleiades a = 0.6 + 0.11 [23], and in cr Orionis 
a = 0.6 + 0.1 ED). 



3.2 Formation radius and the brown dwarf desert 



We have performed an ensemble of 12 radiative hydro- 
dynamical simulation of unstable discs. These discs frag- 
ment on a dynamical timescale (a few thousand years) and 
~ 100 objects are produced overall; mainly BDs but also 
hydrogen-burning stars and planetary-mass objects. 



The discs fragment at distances from 100 to 300 AU from 
the central star. Most of the fragments form with an ini- 
tial mass as low as 3 Mj ( 11321 ). After a fragment forms 
it accretes mass from the disc, and interacts with the disc 
through drag forces and dynamically with other fragments. 
As a result some fragments migrate close to the central star. 
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Fig. 3. The BD desert. The orbital semi-major axes of the objects 
formed by disc fragmentation and remain bound to the central 
star plotted against the mass of the objects. The bars indicate the 
minimum and maximum extent of the orbit. There is a lack of BD 
companions close to central star, but there is a population of BDs 
loosely bound to the central star (at distances ~ 100 - 5000 AU). 

This region is rich in gas, hence these fragments accrete 
material from the disc and eventually become low-mass 
hydrogen-burning stars. Fragments that form farther out 
also accrete material from the disc but not quite as much; 
hence, they remain in the sub-stellar mass regime, i.e. they 
become BDs. If any of the BDs happen to migrate in the 
region close to the central star they tend to be ejected back 
into the outer region through 3 -body interactions. Hence, 
the disc fragmentation model produces a lack of BD close 
companions to Sun-like stars, i.e the BD desert (Fig. [3). 
The BD desert may extend out to 300 AU but it is less dry 
outside 100 AU. Most of the BD companions formed by 
disc fragmentation have eccentric orbits, hence they could 
spend some time within 100 AU from the central star dur- 
ing their periastron passage. 

3.3 Where to look for brown dwarfs forming in 
discs? 

A large fraction of the BDs forming in discs are ejected 
into the field, with typical velocities of ~ 1-3 kms 4 . 
The BDs that remain bound to the central star are on wide 
orbits, i.e. they are loosely bound. Hence, we predict a 
population of BDs at distances from 20 to 5000 AU from 
the central star. This is the region that should be targeted 
by surveys looking for BD companions to Sun-like stars. 
However, in a cluster environment such loosely bound com- 
panions could be liberated (|13|) and become field BDs; 
hence, the number of BDs on wide orbits should be smaller 
than what our simulations of isolated discs predict. 

3.4 Brown dwarf discs 

Most of the BDs form with discs with masses of a few tens 
of Mj and radii of a few tens of AU. These are consistent 
with observations of BD discs (031, lTl9lO . 



3.5 Free-floating planetary-mass objects 

The disc fragmentation model provides an explanation of 
the existence of free-floating planetary-mass objects (|20|, 
[36 1). In our model these objects form in the disc by grav- 
itational fragmentation and quickly after their formation 
they are ejected from the system due to 3-body interac- 
tions. Hence, they stop accreting and their mass remains 
low. Based on the specific initial conditions we examine, 
we predict that free-floating BDs outnumber free-floating 
planetary-mass objects by a factor of at least 10. 

3.6 Low-mass binary statistics 

Close and wide BD-BD binaries are a common outcome 
of disc fragmentation. The components of the binary form 
independently in the disc and then pair up. The simulations 
produce all kinds of low-mass binaries: star-star, star-BD , 
BD-BD, and BD-planetary mass binaries. 

Binary fraction. Our model predicts a low-mass binary 
fraction of 16%. This is comparable with the low-mass bi- 
nary fraction observed in star-forming regions (e.g. in Tau- 
rus : > 20% 07), in the field 15+5%, 03). However, these 
are optical surveys and they are not able to detect tight 
binaries with separation of a few AU. Very tight binaries 
can be probed by radial velocity surveys (e.g. ifToTl ). which 
find a low-mass tight-binary fraction of 10 - 30%. This 
is roughly consistent with our model. However, our model 
produces either very tight or very wide binaries (|30|). 

Wide and close low-mass binaries. We predict that 
close low-mass binaries should outnumber wide ones, and 
this seems to be what is observed (10). A few of the low- 
mass binaries, including one of the very wide systems, 
have become unbound from the central star. Therefore the 
ejection mechanism does not militate against delivering 
wide low-mass binaries to the field. Amongst the low-mass 
binaries that remain bound to the central star, the semi- 
major axis a of the orbit around the central star tends to in- 
crease as the total mass of the low-mass binary decreases. 

Mass-ratios of low-mass binaries. Most of the low- 
mass binaries (55%) have components with similar masses 
(q > 0.7), in agreement with the observed properties of 
low-mass binaries (e.g. 0). 

Low-mass binaries as companions to Sun-like stars. 

The model also predicts that BDs that are companions to 
Sun-like stars are more likely to be in binaries (binary fre- 
quency 25%) than BDs in the field (frequency 5 to 8%). 
This trend is comparable to what is observed, although the 
observed binary frequencies are somewhat higher. [8| re- 
port a binary fraction of 45*}?% for BD companions to 
Sun-like stars, and a binary fraction of only 18^% for BDs 
in the field. This is corroborated by the results of [ 10 1. 
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4 Can planetary-mass objects form by disc 
fragmentation? 

A large number of exoplanets has been observed around 
Sun-like stars. It has been suggested that these planets could 
form by gravitational fragmentation of discs (J7), J5|). In 
our simulations most of the fragments initially have masses 
as low as 3 Mj ( 113211 . see also O). These fragments form 
in the outer disc region (> 100 AU) as the inner disc re- 
gion cannot cool fast enough ( l29l ). However, these ob- 
jects quickly grow in mass by accreting material from the 
disc. Only fragments that are quickly ejected from the disc 
after their formation end up with masses below 13 Mj. 
Hence, it is unlikely that the exoplanets observed close to 
Sun-like stars form by gravitational fragmentation of discs. 



5 Observing fragmenting discs 

Typically a few BDs form in each fragmenting disc around 
a Sun-like star. Hence, it may be that only ~ 20% of Sun- 
like stars need to have large and massive unstable discs to 
produce a large fraction of the observed BDs. Assuming 
that the lifetime of the Class phase is 10 5 yr and that 
in the disc fragmentation scenario the disc fragments and 
therefore dissipates within 10 4 yr , then the probability of 
observing a fragmenting disc around a Class object is 
only 10%. Then, considering that only 20% of Sun-like 
stars may have such unstable discs, the probability of ob- 
serving such discs is only 2%. Hence, fragmenting discs 
should be very difficult to find (||22|), but there are hints of 
their existence (HI)- 



6 Conclusion 

The mechanism of BD formation by fragmentation of grav- 
itationally unstable discs reproduces the BD desert around 
Sun-like stars, and can explain the existence of free-floating 
planetary mass objects and the binary properties of low- 
mass objects. We suggest that there is a population of BDs 
on wide orbits around Sun-like stars that could be probed 
with future observations. 
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